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I. INTRODUCTION 



Progress in the field of ultracold molecules has been rapidly growing for the past decade, 
as many atomic and molecular research groups have turned towards the study of either 
production or dynamics of ultracold diatomic molecules. The success behind the study of 
ultracold molecular formation lies in the use of photoassociation and Feshbach resonances 
(see Jones et al.-^ and Kohler et al.^ for recent reviews). Using a combination of photoas- 
sociation and STIRAP^ (stimulated Raman adiabatic passage) the formation of vibrational 
ground state KRb^i^ and Cs2r molecules have been reported. With this prospect of ground 
vibrational diatoms in mind we focus our study to the interaction effects of a colliding 
lithium atom with that of a f = singlet lithium diatom. Knowledge of both the long- and 
short-range^ interaction potential can be used to calculate inelastic three body rates* as well 
as investigate Efimov coUisions.- 

The assumption of a vibrational ground state diatom greatly simplifies the physics of the 
interaction, however the complexity of calculating long-range interactions for even the f = 
range of motion of the diatom is significant. Given this we work within the frozen diatom 
approximation, which entails freezing the diatom at the calculated equilibrium bond length 
for the entire potential energy surface, which is analogous to vibrationally averaging the 
diatom over the course of the collision. The long-range interaction is not strongly effected 
by this approximation due to the small linear changes seen in the diatomic polarizability 
as the diatom undergoes small oscillations. Here then the polarizability averaged over the 
vibrational motion is just the polarizability evaluated at r^. The validity of the rigid rotor 
approximation for short atom-diatom coUisional distances was found to be good for distances 
greater than lOA. 

We structure this work into the following three parts, first is the discussion of the ah initio 
calculations. Here we discuss the methodology involved in choosing the appropriate basis set 
which provides both the optimal values for the spectroscopic constants as well as the best 
atomic static polarizability. Further improvement to the interaction energy is shown to come 
from accounting for basis set superposition error through a counterpoise correction^- and 
from extrapolating the counterpoise corrected energies to the complete basis limit. Next we 
present a summary of the long-range van der Waals interaction for the tri-atomic dissociation 
limit. An account of the many-body terms arising from third and fourth order Rayleigh- 
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Schrodinger perturbation theory that contribute to the three-atom Cq and Cg van der Waals 
coefficients is given. Finally we present the long-range analytic van der Waals atom-diatom 
interaction energy of Cvitas et ai— and our fitted non-additive van der Waals coefficients. 

II. AB INITIO CALCULATION 

The long-range potential energy surface for the "^A' Lis state has been calculated within 
a frozen diatom approximation for collisional angles near the C2v geometry, which corre- 
sponds to a lithium diatom in the state colliding with a single ^5* state lithium atom. 
The long-range potential energy surface was calculated for atom-diatom distances ranging 
from lOA to lOOA, a region we consider to be both outside of any consideration of charge 
overlap yet still well within the very-long-range limit where the retarded potential starts con- 
tributing. The collisional angle sampling of 6' = 60, 70, 80 and 90 degrees (where 6 = 90° 
corresponds to C2v geometry in Jacobi R, r, 6 coordinates) is both sufficient for an accurate 
fit and consistent with our previous work on the near equilibrium geometry potential energy 
surface.- All electronic energy calculations in this work were done correlating all electrons 
using spin restricted coupled cluster theory with singles, doubles and iterative triples^^*^ 
(RCCSD(T)) as implemented in the MOLPRO 2008.1— suite of ab initio programs. 

Cold scattering calculations require exceptionally accurate interaction potentials in order 
to properly predict cross-sections and scattering lengths. To provide this accuracy we apply a 
series of corrections to the RCCSD(T) energy which account for deficiencies within the basis 
set. Additionally we correlate all electrons in the RCCSD(T) energy so as to account for 
the core-core and core-valence (CV) contributions. The inclusion of CV correlation energy 
has been shown to account for roughly 0.002A for multiple bonds and several hundred cm~^ 
to atomization energies.— To properly correlate all electrons within an ab initio calculation 
necessitates the use of a CV consistent basis set.— We have examined the use of the four 
and five zeta cc-pVnZ basis sets from Feller— and the CV consistent CVnZ basis sets of 
Iron et al}^ In Table H] is a comparison between calculated spectroscopic constants from the 
above basis sets and the experimental values, it can be seen that the five zeta basis sets 
provide a marked improvement in terms of both rg and Dg. The calculated polarizability is 
also used as a benchmark calculation in addition to and Df,. The accuracy of which is a 
strong indicator of the accuracy of the atomic static polarizability is a strong indicator of 
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the accuracy of long-range molecular dispersion interactions. 

The first improvement to the trimer potential energy surface was to calculate the basis 
set superposition error (BSSE) through a counterpoise calculation.— The interaction energy 
of the trimer is then 

rp rp pis pis P^S /-i \ 

-C/int — -C/trimer -'^atom -'^atom -^atom 

where E]^^^^ is the energy of the resulting atom when atoms i and j within the trimer are 
replaced with dummy centers. Accounting for BSSE the CVQZ basis provides a correction 
of 5.27 cm^^ and the CV5Z basis has a correction of 1.05 cm~^. Across the total long-range 
surface the total BSSE correction varies no more than 1%, finally converging to a constant 
value at interaction regions greater than 20A. This suggests that the majority of the BSSE 
corresponds to the diatom and not the atom-diatom interaction itself. Still this counterpoise 
correction does not fully account for all of the dissociation energy of the Li2 diatom. We 
further improve upon the accuracy of the potential energy surface by extrapolating the 
counterpoise corrected energies to the complete basis set (CBS) limit. We use the CBS limit 
extrapolation formulation of Helgaker et ai— 

En = Eoo H — ^. (2) 

This extrapolation scheme was applied to the CVQZ and CV5Z counterpoise corrected 
interaction energies for the final potential energy surface calculation. 

We calculate the static polarizability for the ground state lithium atom and singlet lithium 
diatom with the static field methocU^iSS by calculating the RCCSD(T) energies in the pres- 
ence of an electric dipole field (given here as E{F)). The static polarizability is given by the 
finite field gradient 

5E{0)/2~8E{F)/3 + E{2F)/6 
a^^ = (3) 

and reported in Table [TTl The dispersion energy between two monomers at long-range (no 
charge overlap) can be expressed as^ 

^ 3Vaia2 , 
Udisp — , (.4j 

where is the static polarizability for the ith monomer and V is a. characteristic excitation 
energy of the molecule^ . Due to the important contribution of the long-range tail of the 
electron wave function to the molecular polarizability, the effects of adding diffuse functions 



can be significant. In the calculation of the polarizability discussed above, a set of even 
tempered diffuse functions were added to the CVnZ basis sets and found to contribute little 
to the extrapolated polarizability. A further test on the effect of adding diffuse functions 
was performed by evaluating a CBS extrapolated single point energy calculation at lOA with 
the even tempered diffuse functions discussed above. The results showed that in the CBS 
limit the difference between the standard and augmented CVnZ basis sets amounts to less 
than half a wavenumber. 

Because the long-range interaction depends explicitly upon the monomer static polariz- 
ability, it is clear from the reported atomic static polarizability in Table [III that the CVnZ 
basis sets are the optimal. The variation in theoretically reported parallel and perpendicular 
polarizabilities for the lithium diatom (see Deiglmayr et af^ for a review of alkali diatomic 
static polarizabilities) do not offer a significant constraint on the chosen basis set. Our final 
recommended calculations make use of the counterpoise correction with the extrapolation 
to the CBS limit as discussed above. We expect that our calculated diatomic long-range in- 
teractions will be suitably accurate given the precision of the calculated spectroscopic values 
and static atomic polarizability. 



III. LONG-RANGE VAN DER WAALS POTENTIAL 

We now examine the trimer potential energy surface at the three atom dissociation limit, 
which is expanded analytically in terms of the two-body dispersion interactions with the 
addition of a purely three-body interaction potential. The three-body contribution to the 
interaction energy is well known to be strong for Lis, for both the doubleti^ and quartet^ 
state. Thus it is important to accurately include such effects in any model of the long- 
range interaction. By using perturbation theory to examine this expansion it is possible 
to express the generalized dispersion coefficients in terms of known diatom and triatomic 
constants. Our goal in this section is to overview the essential theory of three-body atomic 
interactions, which can then be specialized to the case of atom-diatom interactions within 
the frozen diatom approximation previously discussed. 

The tri-atomic dissociation interaction potential can be described in terms of the diatomic 
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van der Waals interaction potential Vd{rij) and non-additive many-body potential V3(r) as 

Vir) = J2v,ir,,) + Vsiv), (5) 

where r = (ri2, ri3, r23) are the three internuclear vectors. Long-range dispersion interaction 
potentials (excluding retardation and orbital overlap effects) between two S-state atoms are 
described using the multipole expansion 

V,{r,,) = -C,rr^ - C,rr^ + O(rriO), (6) 

where the Cq and Cs coefficients are respectively the dipole-dipole, quadrupole-quadrupole 
and dipole-octopole expansions^i of the inter-atomic electrostatic Hamiltonian Vij. To ob- 
tain the leading terms of the non- additive potential V3(r) in Eq. El the analogous expan- 
sion method used for the diatom van der Waals interaction can be used. Here, Rayleigh- 
Schrodinger (RS) perturbation theory is applied to the total inter-atomic interaction Hamil- 
tonian ifint = H^o^ — Ha — Hb — Hci where is the atomic Hamiltonian for the x'th atomic, 
then expanded by multipole moments. The desired contributions to the non-additive poten- 
tial V3(r) arise from the first many-body terms in third order RS perturbation theory. 

The bipolar expansion of the many-body terms from third order RS perturbation theory 
leads to a summation of purely geometric factors with interaction constants dependent on 
the atomic speciesS^iS^, 

V^'^ = Y.^^mMZim,. (7) 

The interaction constant can be expressed with the Casimir-Polder integral of the dynamic 
2'' polarizabilities^^ over complex frequencies 

1 /"°° 

^hhh = ~ a^^{iuj)a^^{iuj)a^'^{iuj)duj. (8) 
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With the dynamic polarizability defined as^ 



n 

and the 2' oscillator strength defined as 



/S = ^S„|{0| 5]r;y;„,(r-,)|n)|^ (10) 



the sum over i goes over all the electrons in the given atom and En is again the excitation 
energy for state n. The geometric factors Wi-^i^i^lr) have been reported by a number of 
authorsiii^^"— and will not be reproduced here. The first term in the third order expansion 
can be identified as the well known Axilrod-Teller-Muto triple-dipole^ term. Additionally, 
it has been identified that the quadruple dipole term Zmi, from fourth order perturbation 
theory, has a contribution to the van der Waals dispersion coefficients of consideration here. 
This term has no exact expression in terms of the dynamic polarizabilitiesi^ , but it can be 
approximated using a Drude oscillator model in the case of three S state atoms. Using the 
corresponding Drude oscillator approximation for the Cg van der Waals dispersion coefficient, 
the Ziiii term can be approximated as^S 

Znii = 10ZIJ{3C,). (11) 

Using the values Cq = 1393.39 and Zm = Vabd'i = 56865 from Yan et al.—, this provides 
the approximate value Zuu = 7735638 a.u. for the quadrupole dipole term. 

In the other asymptotic limit of the trimer where the system dissociates to a diatom and 
separated atom, the van der Waals type interaction can again be expressed as a series of 
multipole terms of the diatomic and atomic polarizabilities. Using the Jacobi coordinates 
R, r and 6 to describe the atom-diatom system, where r is the diatomic internuclear distance, 
R is the diatomic center of mass to colliding atom distance and 6 is the angle between r and 
R, for the asymptotic limit of R ^ r the interaction potential in the absence of damping 
and exchange is^i 

V{R, r, 9) = D,- Ce{r, e)R-^ - Cs{r, e)R-'' + 0{R-'') (12) 

Where is the interaction energy of the diatom, and the dispersion coefficients are defined 
in terms of Legendre polynomials as 

C,ir,e) =C°(r) + Clir)P,{cose), (13) 
Cs{r,9) =C's{r) + Cl{r)P2{cos9) + Cl{r)P,{cos9) (14) 

(15) 

To obtain the analytic form for the atom-diatom van der Waals coefficients in terms of 
the tri-atomic terms, we transform the internuclear rij coordinates to Jacobi coordinates 

7 



through the transformations 

ri2 = r, (16) 

2 

r23 = (R'' + J + Rr COS ey/^ (17) 

2 

rgi = (i?2 + L _ COS 0)^/^ (18) 

in the hmit of r. 

The contributions to Cq^v, 0) and C^{r, 6) can be found by transforming the Wi^i^i^ir) and 
coefficients and then expanding in a power series with respect to r/R. This expansion 
has been completed by Cvitas et al}^ for the contributions to C^ir^O) and C^{r^O). The 
resulting terms included in these van der Waals coefficients are 

C°(r) =2C6 + 12Zniir-6 + ©(r'^), (19) 
Cl{r) =6Znir-3 + GZniir-^ + 0(r-8), (20) 

C°(r) =2C8 + ^^r^ - ^Znir-'+ 

yZnnr-^ + 0(r-6), (21) 

^°^-^iiiir-^ + 0(r-6) (22) 
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and 



(40Zn3 - 30Zi22)r-^ + 0(r-6). (23) 

In the atom-diatom dissociation limit the effects of charge overlap damping in the r 
dependent terms of Eqs. [T9ll23l cannot be neglected. The standard method of accounting 
for charge overlap is the inclusion of a damping coefficient fn{fij) for each r^^ contribution. 
Recent reports^ii^ have made use of the popular Tang and Toennies^ dispersion damping 
functions 

/n(r,) = l-e-^-^X:^- (24) 

fc=0 

In this work we chose to follow Rerat and Bussery-Honvault22,'s implementation of the Tang 
and Toennies damping functions, where each anisotropic contribution in terms of r has the 
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following mapping 

r-" ^ /n(r)r-". (25) 

Working within the frozen diatom approximation the charge overlap damping can be simply 
modeled as a constant fitting parameter, evaluated at the diatomic equilibrium bond 
length. Inserting the fitting parameter into Eqs. [T91I23] provides the following 

C°(re) =2C, + l2Z^^^^F,rf + 0(r-«), (26) 
Cl{r,) =6ZniFsr:^ + QZ^mF^rf + ©(r'^), (27) 

C0(r.)^2C3 + ^CeF...^-^W.r-+ 

^Z,,,,F,T:' + 0{r~% (28) 
C|(r,) =SC,F2rl + ^Z^iiF^r'^ + l|^Zn2F3r;3+ 

^^niiF4r,-^ + 0(r-^) (29) 

and 

Ct[r:) =^-^ZniF,r;' - ^^ii2F3r,-3 + ^-^ZnnF,r:'+ 

(40^113 - 30Zi22)F5r;' + 0(r-6). (30) 

As evaluated in Eqs. I26ll30| the fitting parameter F„ no longer directly correspond to the 
Tang and Toennies damping function as has been noted.— 

We have performed a linear least squares fit of Eq. [121 with the definitions given by Eqs. 
[T3]fT^ to obtain the van der Waals coefficients given in Table IIIII We have also fitted Eq. [12] 
with the definitions given in Eqs. [TU1I2^ We use the diatom Cg = 1393.39 and Cg = 8342 
dispersion coefficients as well as the Zm = Vabd^ = 56865 triple dipole term from Yan et 
alM and the Zu2 = 581000, Zu3 = 1.7 x 10^ and Z122 = 6.41 x 10^ from Patil and Tang.^ 
Fitting to these values we obtain the same dispersion coefficients as above due to the nature 
of the least squares fit, with the values for the damping parameters given in Table [XTTl Plotted 
in Fig. [T] is a comparison of the ab initio surface and the fitted van der Waals potential. 
As can be seen the fit is very accurate, which is confirmed by the RMS surface fitting error 
of 10~^ cm~^. This analytical long-range expansion can be used in conjunction with our 
previous work'' on the ground state surface of "^A' Lis to evaluate scattering properties of 
the Li+Li2 rigid rotor system. The results of which we leave to future work. 



IV. CONCLUSIONS 



We have calculated ab initio an accurate long-range "^A' Lis surface for the dissociation 
to the Li pS]+ Lia [X'^T.+]. The lithium diatom was taken to be rigid rotor with the bond 
length constrained to r = re, the calculated equilibrium bond length in table HI The surface 
was calculated at the RCCSD(T) level of theory, correlating all electrons. To accurately 
describe the CV interaction, the CVQZ and CV5Z basis sets of Iron et al}^ were used; the 
final counterpoise corrected interaction energies were then extrapolated to the CBS limit. 
At this level of theory, the atomic and diatomic dipole polarizabilities are shown to be in 
good agreement with published experimental and theoretical results, which is an important 
component to long-range interactions. Using the expansion for the atom-diatom many body 
van der Waals potential (Eq. [T2]and [T3l) we calculate the non- additive interaction coefficients 
by fitting to the calculated ah initio surface. The resulting repulsive contribution of the 
three-body interaction to the total interaction energy is found to be up to 33% of the total 
energy. The fitted van der Waals coefficients were found to be consistent with the existing 
literature on the related '^A' state.— 
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basis 


A Dq cm ^ 


Df. cm ^ 


Exp. 


2.6731 


8434.58^. 


8516.78^ 


Extrapolated 


2.674 




8510.926 


CV5Z 


2.674 


8310.85 


8487.402 


CVQZ 


2.676 


8288.66 


8465.005 


cc-pV5Z 


2.680 


8319.12 


8470.085 


cc-pVqZ 


2.685 


8270.86 


8445.741 



^ reference^ 
^ reference^ 

TABLE I. Spectroscopic constants for the ground state singlet lithium diatom calculated using 
four and five zeta correlation consistent basis sets with counterpoise. 



Li pS] Li2 



basis 



a an 



a_L aj] an aj] 



Extrapolated 164.1 
CV5Z 164.2 
CVQZ 
cc-pV5Z 
cc-pVqZ 
Theory 



164.3 
165.6 
166.7 
164. IIP 
164.4^ 
Experiment 164.2(1.0)^ 



165.6 296.1 

165.1 296.7 

164.6 297.2 

166.4 299.6 

166.4 300.8 

162.4^ 305.2^ 



reference^i 
^ reference^^ 
reference^ 

TABLE II. Calculated finite field static polarizabilities for the ground state lithium atom and 
singlet lithium diatom evaluated at the RCCSD(T)/CV5Z equilibrium bond length. 
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D ^ 


= 1604 


Fi = 


-0.2823 




= 120 


F2 = 


0.1602 




= 221646 


F3 = 


0.0399 


Ciire) 


= 173870 


Fa = 


0.0167 


Ciire) 


= 27104 


F5 = 


0.0075 






F6 = 


-0.0047 



TABLE III. Fitted van der Waals coefficients in a.u. for the Li-|-Li2 ridged rotor system. 




FIG. 1. Comparison of the ab initio surface (points) and fitted van der Waals potential (solid lines) 
for the Li [^S]+ Li2 [X'^T,^] interaction. The Li2 bond length is held at the calculated equilibrium 
bond length rg. 
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